Electronic Correlations in Double Quantum Dots 
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We present a study of the electronic structure of two laterally coupled Gaussian quantum dots 
filled with two particles. The exact diagonalization method has been used in order to inspect 
the spatial correlations and examine the particular spin singlet-triplet configurations for different 
coupling degrees between quantum dots. The outcome of our research shows this structure to have 
highly modifiable properties promoting it as an interesting quantum device, showing the possible 
use of this states as a quantum bit gate. 
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The growing significance of mesoscopic quantum dots 
(QDs) in the development of microelectronic devices at- 
tracts great attention, partly due to their promising tech- 
nological utilization, but also due to their novel quantum 
properties. These two-dimensional nanostructures offer 
a high degree of flexibility in their manufacturing pro- 
cess giving as a result a very interesting combination of 
a highly controllable structure showing quantum effects. 
Recent advances in the manipulation of electrically de- 
fined QDs open a new line of investigation in coupled 
quantum dots (CQDs) and a very interesting chance for 
the realization of a solid state implementation of the ba- 
sic components for quantum computing. The interest on 
CQDs in the field of quantum computing was brought 
over since Loss and DiVincenzo proposal^ on the pos- 
sibility of an implementation based on the spin states 
of coupled single-electron quantum dots. Recent studies 
and experiments^i^i^i^i^ have proven the feasibility of us- 
ing the spin degree of freedom in vertically and laterally 
coupled few-electron QDs as systems in which interac- 
tions can be electro-optically tuned and therefore they 
can potentially taken as candidates for performing quan- 
tum operations. 

It is the purpose of the present paper to present an ac- 
curate theoretical study focused on the electronic struc- 
ture of the two-particle GDQD in the absence of a mag- 
netic field. The exact diagonalization method^ has been 
used in order to solve the interacting many-body system 
as it is a tool of proven reliability for few-particle prob- 
lems^. Electron-electron correlation effects are studied 
for the different spin superposition states and for a vari- 
ety of configurations, from a weakly interacting laterally 
coupled quantum dot molecule to a strongly overlapped 
situation. The choice of an anisotropic Gaussian confin- 
ing potential makes it possible to describe a more real- 
istic situation.* We will show that not only in vertically 
coupled QDs interactions can be tuned but also in lat- 
eral QDs this effect can be achieved; furthermore, this 
systems could actually be used as dynamically tunable 
two-level systems. 

As shown in our previous research^, the special be- 
havior of a GDQD filled with two electrons makes it a 



very promising arrangement for the creation of a dynam- 
ically tunable quantum gate. As noted by previous re- 
search*-^'^'^, one of the most attractive setup among few- 
particle lateral quantum dots is the two-electron case in 
which the singlet and triplet states can be tuned to be 
bound into a fourfold degenerate ground state or, de- 
pending in the double QD parameters, artificially intro- 
ducing an energy gap between the singlet and triplet 
states. This particular characteristic is very similar to 
that already observed in a two-minima QD subjected to 
a magnetic field^. The possibility of using two coupled 
quantum dots containing two electrons as a charge qubit 
by dynamically tuning the interdot interaction has been 
recently proposed by Weiss et ali^. Despite them using 
two vertically coupled QDs, our system is physically sim- 
ilar and it could also be employed in developing a charge 
qubit. 

In our model, the single-particle two-dimensional 
Hamiltonian for a particle subjected to a magnetic field 
is given by 
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where m* is the effective mass of the electron, e = — |e| is 
the electron charge, g* is the effective Lande g-factor and 
VG(r) is the external confining potential. We have used 
the model by Hu and Das Sarma* in order to describe 
a realistic trapping potential. In this case, Vg is formed 
by three anisotropic Gaussians, two of them modeling 
the lateral quantum dots and the third one controlling 
the height of the potential barrier between them. More 
precisely, VQ{x,y) can be written as 
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with all parameters being positive except Vq, which con- 
trols the strength of the lateral QDs. xq regulates their 




FIG. 1: Electron density contours for a GDQD with Vb = 4.0 meV (left), 8.5 meV (center) and 14.0 meV (right). Left and 
center configurations show a clear coupling between both dots, whereas circular symmetry in each one of the dots is practically 
restored when Vb is set to 14.0 meV. The point of maximum density is always placed at y = due to symmetry and, along the 
X-axis is placed at ±24 nm, ±37 nm and ±45 nm from left to right. These points will be used when calculating conditional 
probabilities. We must note that setting Vb to a sufficiently large and negative value, so that Vb is negligible compared with 
Vb, we could restore a case in which only a single minima Gaussian QD is obtained. 



relative distance. As we have mentioned before, we must 
note that the QDs studied in this paper do not exhibit 
circular symmetry as, in general, ^ ly and Ihx hy In 
this model, the parameter VJ, independently modulates 
the height of the interdot potential barrier, effectively 
controlling the interdot distance. It is important to note 
that solely varying the parameter V;, we can smoothly 
jump from a single quantum dot filled with two electrons 
or, in the opposite case, with two independent quantum 
dots each of them occupied by a single electron. This 
is a very interesting parameter as it could be dynam- 
ically modified in an electrically defined QD giving the 
chance to exploit all of its properties. For the many-body 
Hamiltonian, it will be enough to include the Coulomb 
electrostatic potential, given by 

^c(ri,r2) = ^r-^ (3) 

We have used the exact diagonalization (ED) method 
in order to solve the Hamiltonian problem. Under this 
approach, we can write the eigenfunctions in the form of 
two-component spinors 

ip = ^c^,,s(t)f,x', (4) 

K,S 

where and ^I'e the eigenvectors of the Sz operator 
and the index k points to a set of quantum numbers. The 
spatial wave functions chosen are the complete set of the 
solution of a two-dimensional anisotropic harmonic oscil- 
lator with characteristic lengths Ux and ay. The many- 
body interacting Hamiltonian is formed by including the 
two-body matrix elements of the Coulomb interaction, 
in the same fashion as in Refj^. Given this basis, the 



Hamiltonian matrix elements can be written in a closed 
and analytic form.^ 

Using the second quantization formalismi^ one can de- 
fine the two-body density operator or the pair correlation 
function as 

p(xi,X2) = ^(/>*(xi) (/.*(x2)(/>fe(x2 )(?!); (xi) (5) 
ijkl 

where x = {r, (} refers simultaneously to spatial and spin 
variables. If we are interested only in spatial correlations, 
we can sum over the spin degree of freedom to obtain 

p(ri,r2) = pn(ri,r2) ±PTi(i"i,i"2) 

+Plt(ri,r2) + /3xi(ri,r2), (6) 

where p||(ri,r2) is the pair correlation function for for 
two spin up electrons and /9||(ri, r2) is the same quantity 
but for two anti-parallel particles. 

It is our purpose to investigate the electronic struc- 
ture of this system under different conditions. For ac- 
complishing this objective, we have chosen the follow- 
ing set of shared parameters for the confining potential 
Vg- Vo — 50 meV, xq = 54 nm, = ly = 80 nm, 
hx — 40 nm and Iby — 50 nm. Material parameters are 
those appropriate of GaAs, m* = 0.067mo, e' = 13.1, 
g* = —0.44. Given all these values, we have done dif- 
ferent calculations for the parameter VJ,, from a strongly 
bound double QD in which VJ, — 4.0 meV to a case in 
which both QDs are nearly isolated from each other be- 
ing Vi, — 14.0 meV, going trough an intermediate con- 
figuration in which Vf, = 8.5 meV. We will consider this 
intermediate setup as our base case. For this case, we can 
see in Fig. [5] the energy dependence on magnetic field for 
a single- and double-electron configuration. There, we 
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FIG. 2: Energy dependence with magnetic field for a single- 
(left) and two- electron (right) GDQD. Interdot potential 
strength is fixed to Vt = 8.5 meV. The two lowest energy 
branches in the case of a single-electron are formed by two 
levels each, merging as magnetic field increases. 

can appreciate that in the single-particle case, ground 
state at B = T is a doubly degenerate level due to spin 
degree of freedom. This degeneracy is temporarily bro- 
ken by magnetic field until reaching a sufficiently large 
B, which couples ground state and second degenerate 
state and makes them degenerate again. In the case of a 
two-electron GDQD, we can see that magnetic field splits 
level spacing monotonically. It is difficult to appreciate 
by naked eye, but energy difference between the first four 
states at i? = T is of the order of the numerical pre- 
cision used for this simulation. If we consider the same 
quantity being Vb = 14.0 meV, this difference will be well 
below our numerical precision, but for the case Vb — 4.0, 
the energy difference is not inappreciable anymore. 

The total electron density p(r) — px(r) + p|(r) for the 
ground state is plotted in Fig. [T] We can appreciate a 
clear separation of the quantum dots as Vb is increased, 
in fact, in the case for higher Vb, we almost recover two 
clearly differentiated and circularly symmetric QDs. In- 
terestingly, the effect in electron density when modifying 
the parameter Vb is similar to that observed by a vari- 
able magnetic field noted by Mario and Helle^'^"^. From 
these figures we can find out the maximum electron den- 
sity points in order to use them when calculating the 
correlations between electrons. Due to axial symmetry, 
these point will be placed at y = nm, and they will be 
symmetric with respect to x-axis. From left to right, the 
abscissa is 24, 37 and 45 nm. Given this values, we have 
plotted Figs. H and H 



Results for electron-electron correlations when Vb = 
8.5 meV are plotted in Fig. [31 One can appreciate a 
clear symmetric-antisymmetric pattern in the distribu- 
tion of the energy levels, this can be explained by the 
spin change of parity in the spinors. For example, in 
Fig. [3] (a) , we can find the conditional probability for the 
ground state, it shows a spatially symmetric distribution, 
this is due to the equal contribution of pn and terms 
into the total correlation /9(ri,r2) while the numerical 
contribution of the remaining terms is numerically negli- 
gible. This symmetry is also observed in subfigureO (c), 
corresponding to the second excited state. In opposition. 
Figs. [3] (b) and [3] (d), those corresponding to first and 
third excited states, are fundamentally composed by the 
spin-symmetrized and pn terms. As we are deal- 
ing with an even ground state has a spatially symmetric 
distribution as well as second excited state. 

The energy spread between the first four lowest ly- 
ing levels, i.e. singlet and triplet states, varies dra- 
matically as a result of modifying the interdot barrier 
height.^ Beginning with the case Vb = 4.0 meV in which 
AE = 65.30 peV and ending at Vb = 14.0 meV we find 
that AE = 0.12 peV. In the intermediate step, Vb — 8.5 
meV the energy spread is 5 peV. The numerical error 
derived from our diagonalization routines for the highest 
barrier case is estimated to be of the same magnitude as 
the energy difference between singlet and triplet states, 
this impedes us to clearly discern between the four en- 
ergy levels and indicates us that we must take them as a 
group of degenerate levels. This situation contrasts with 
the cases with a lower interdot barrier, in which energy 
separation is clear and there is no evidence of energy level 
degeneracies to appear. 

In summary, we have studied the effect of electrostatic 
interactions in a two-electron laterally coupled double 
quantum dot. The results of our exact diagonalization 
calculations show a peculiar singlet-triplet level merging 
as coupling strength is screened by a interdot barrier. 
Further analysis of the electronic wave functions such as 
calculating conditional probabilities of the lowest energy 
levels, show the change produced in the spatial and spin 
degrees of freedom. 
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FIG. 3: Pair-correlation function p(r, r') for the GDQD with Vt — 8.5 meV and r' — {x',y') = (37,0) nm. Figure (a) shows 
ground state conditional density and figures (6), (c) and (d) correspond to first, second and third excited states. Scale in x and 
y axes is measured in nm. 




FIG. 4: Pair-correlation function p(r, r') for the GDQD with Vb = 14.0 meV and r' — {x', y') — (45, 0) nm. Figure (a) shows 
ground state conditional density and figures (6), (c) and (d) correspond to first, second and third excited states. Scale in x and 
y axes is measured in nm. 
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